Transcriptional and metabolic profiling of sulfur starvation response in two monocots.

BACKGROUND: Sulfur (S) is a mineral nutrient essential for plant growth and development, which is incorporated into diverse molecules fundamental for primary and secondary metabolism, plant defense, signaling, and maintaining cellular homeostasis. Although, S starvation response is well documented in the dicot model Arabidopsis thaliana, it is not clear if the same transcriptional networks control the response also in the monocots. RESULTS: We performed series of physiological, expression, and metabolite analyses in two model monocot species, one representing the C3 plants, Oryza sativa cv. kitaake, and second representing the C4 plants, Setaria viridis. Our comprehensive transcriptomic analysis revealed twice as many differentially expressed genes (DEGs) in S. viridis than in O. sativa under S-deficiency, consistent with a greater loss of sulfur and S-containing metabolites under these conditions. Surprisingly, most of the DEGs and enriched gene ontology terms were species-specific, with an intersect of only 58 common DEGs. The transcriptional networks were different in roots and shoots of both species, in particular no genes were down-regulated by S-deficiency in the roots of both species. CONCLUSIONS: Our analysis shows that S-deficiency seems to have different physiological consequences in the two monocot species and their nutrient homeostasis might be under distinct control mechanisms.

Variations and reduction of plastome are associated with the evolution of parasitism in Convolvulaceae.

Parasitic lifestyle can often relax the constraint on the plastome, leading to gene pseudogenization and loss, and resulting in diverse genomic structures and rampant genome degradation. Although several plastomes of parasitic Cuscuta have  been reported, the evolution of parasitism in the family Convolvulaceae which is linked to structural variations and reduction of plastome has not been well investigated. In this study, we assembled and collected 40 plastid genomes belonging to 23 species representing four subgenera of Cuscuta and ten species of autotrophic Convolvulaceae. Our findings revealed nine types of structural variations and six types of inverted repeat (IR) boundary variations in the plastome of Convolvulaceae spp. These structural variations were associated with the shift of parasitic lifestyle, and IR boundary shift, as well as the abundance of long repeats. Overall, the degradation of Cuscuta plastome proceeded gradually, with one clade exhibiting an accelerated degradation rate. We observed five stages of gene loss in Cuscuta, including NAD(P)H complex → PEP complex → Photosynthesis-related → Ribosomal protein subunits → ATP synthase complex. Based on our results, we speculated that the shift of parasitic lifestyle in early divergent time promoted relaxed selection on plastomes, leading to the accumulation of microvariations, which ultimately resulted in the plastome reduction. This study provides new evidence towards a better understanding of plastomic evolution, variation, and reduction in the genus Cuscuta.

Comprehensive deciphering the alternative splicing patterns involved in leaf morphogenesis of Liriodendron chinense.

Alternative splicing (AS), a pivotal post-transcriptional regulatory mechanism, profoundly amplifies diversity and complexity of transcriptome and proteome. Liriodendron chinense (Hemsl.) Sarg., an excellent ornamental tree species renowned for its distinctive leaf shape, which resembles the mandarin jacket. Despite the documented potential genes related to leaf development of L. chinense, the underlying post-transcriptional regulatory mechanisms remain veiled. Here, we conducted a comprehensive analysis of the transcriptome to clarify the genome-wide landscape of the AS pattern and the spectrum of spliced isoforms during leaf developmental stages in L. chinense. Our investigation unveiled 50,259 AS events, involving 10,685 genes (32.9%), with intron retention as the most prevalent events. Notably, the initial stage of leaf development witnessed the detection of 804 differentially AS events affiliated with 548 genes. Although both differentially alternative splicing genes (DASGs) and differentially expressed genes (DEGs) were enriched into morphogenetic related pathways during the transition from fishhook (P2) to lobed (P7) leaves, there was only a modest degree of overlap between DASGs and DEGs. Furthermore, we conducted a comprehensively AS analysis on homologous genes involved in leaf morphogenesis, and most of which are subject to post-transcriptional regulation of AS. Among them, the AINTEGUMENTA-LIKE transcript factor LcAIL5 was characterization in detailed, which experiences skipping exon (SE), and two transcripts displayed disparate expression patterns across multiple stages. Overall, these findings yield a comprehensive understanding of leaf development regulation via AS, offering a novel perspective for further deciphering the mechanism of plant leaf morphogenesis.

Functional differentiation of olive PLP_deC genes: insights into metabolite biosynthesis and genetic improvement at the whole-genome level.

KEY MESSAGE: This study identified 16 pyridoxal phosphate-dependent decarboxylases in olive at the whole-genome level, conducted analyses on their physicochemical properties, evolutionary relationships and characterized their activity. Group II pyridoxal phosphate-dependent decarboxylases (PLP_deC II) mediate the biosynthesis of characteristic olive metabolites, such as oleuropein and hydroxytyrosol. However, there have been no report on the functional differentiation of this gene family at the whole-genome level. This study conducted an exploration of the family members of PLP_deC II at the whole-genome level, identified 16 PLP_deC II genes, and analyzed their gene structure, physicochemical properties, cis-acting elements, phylogenetic evolution, and gene expression patterns. Prokaryotic expression and enzyme activity assays revealed that OeAAD2 and OeAAD4 could catalyze the decarboxylation reaction of tyrosine and dopa, resulting in the formation of their respective amine compounds, but it did not catalyze phenylalanine and tryptophan. Which is an important step in the synthetic pathway of hydroxytyrosol and oleuropein. This finding established the foundational data at the molecular level for studying the functional aspects of the olive PLP_deC II gene family and provided essential gene information for genetic improvement of olive.

A membrane associated tandem kinase from wild emmer wheat confers broad-spectrum resistance to powdery mildew.

Crop wild relatives offer natural variations of disease resistance for crop improvement. Here, we report the isolation of broad-spectrum powdery mildew resistance gene Pm36, originated from wild emmer wheat, that encodes a tandem kinase with a transmembrane domain (WTK7-TM) through the combination of map-based cloning, PacBio SMRT long-read genome sequencing, mutagenesis, and transformation. Mutagenesis assay reveals that the two kinase domains and the transmembrane domain of WTK7-TM are critical for the powdery mildew resistance function. Consistently, in vitro phosphorylation assay shows that two kinase domains are indispensable for the kinase activity of WTK7-TM. Haplotype analysis uncovers that Pm36 is an orphan gene only present in a few wild emmer wheat, indicating its single ancient origin and potential contribution to the current wheat gene pool. Overall, our findings not only provide a powdery mildew resistance gene with great potential in wheat breeding but also sheds light into the mechanism underlying broad-spectrum resistance.

Fine mapping and breeding application of two brown planthopper resistance genes derived from landrace rice.

The Brown planthopper (Nilaparvata lugens Stål; BPH) is known to cause significant damage to rice crops in Asia, and the use of host-resistant varieties is an effective and environmentally friendly approach for controlling BPH. However, genes limited resistance genes that are used in insect-resistant rice breeding programs, and landrace rice varieties are materials resources that carry rich and versatile genes for BPH resistance. Two landrace indica rice accessions, CL45 and CL48, are highly resistant to BPH and show obvious antibiosis against BPH. A novel resistance locus linked to markers 12M16.983 and 12M19.042 was identified, mapped to chromosome 12 in CL45, and designated Bph46. It was finely mapped to an interval of 480 kb and Gene 3 may be the resistance gene. Another resistance locus linked to markers RM26567 and 11MA104 was identified and mapped to chromosome 11 in CL48 and designated qBph11.3 according to the nominating rule. It was finely mapped to an interval of 145 kb, and LOC_Os11g29090 and LOC_Os11g29110 may be the resistance genes. Moreover, two markers, 12M16.983 and 11MA104, were developed for CL45 and CL48, respectively, using marker-assisted selection (MAS) and were confirmed by backcrossing individuals and phenotypic detection. Interestingly, we found that the black glume color is closely linked to the BPH resistance gene in CL48 and can effectively assist in the identification of positive individuals for breeding. Finally, several near-isogenic lines with a 9311 or KW genetic background, as well as pyramid lines with two resistance parents, were developed using MAS and exhibited significantly high resistance against BPHs.

Fine mapping of TFL, a major gene regulating fruit length in snake gourd (Trichosanthes anguina L).

Fruit length is a crucial agronomic trait of snake gourd (Trichosanthes anguina L); however, genes associated with fruit length have not been characterised. In this study, F2 snake gourd populations were generated by crossing the inbred lines, S1 and S2 (fruit lengths: 110 and 20 cm, respectively). Subsequently, bulk segregant analysis, sequencing, and fine-mapping were performed on the F2 population to identify target genes. Our findings suggest that the fruit length of snake gourd is regulated by a major-effect regulatory gene. Mining of genes regulating fruit length in snake gourd to provide a basis for subsequent selection and breeding of new varieties. Genotype-phenotype association analysis was performed on the segregating F2 population comprising 6,000 plants; the results indicate that the target gene is located on Chr4 (61,846,126-61,865,087 bp, 18.9-kb interval), which only carries the annotated candidate gene, Tan0010544 (designated TFL). TFL belongs to the MADS-box family, one of the largest transcription factor families. Sequence analysis revealed a non-synonymous mutation of base C to G at position 202 in the coding sequence of TFL, resulting in the substitution of amino acid Gln to Glu at position 68 in the protein sequence. Subsequently, an InDel marker was developed to aid the marker-assisted selection of TFL. The TFL in the expression parents within the same period was analysed using quantitative real-time PCR; the TFL expression was significantly higher in short fruits than long fruits. Therefore, TFL can be a candidate gene for determining the fruit length in snake gourd. Collectively, these findings improve our understanding of the genetic components associated with fruit length in snake gourds, which could aid the development of enhanced breeding strategies for plant species.

An pair of an atypical NLR encoding genes confer Asian soybean rust resistance in soybean.

Asian soybean rust (ASR), caused by Phakopsora pachyrhizi, is a devastating disease that is present in all major soybean-producing regions. The limited availability of resistant germplasm has resulted in a scarcity of commercial soybean cultivars that are resistant to the disease. To date, only the Chinese soybean landrace SX6907 has demonstrated an immune response to ASR. In this study, we present the isolation and characterization of Rpp6907-7 and Rpp6907-4, a gene pair that confer broad-spectrum resistance to ASR. Rpp6907-7 and Rpp6907-4 encode atypic nucleotide-binding leucine-rich repeat (NLR) proteins that are found to be required for NLR-mediated immunity. Genetic analysis shows that only Rpp6907-7 confers resistance, while Rpp6907-4 regulates Rpp6907-7 signaling activity by acting as a repressor in the absence of recognized effectors. Our work highlights the potential value of using Rpp6907 in developing resistant soybean cultivars.

Application of Integrated Computational Approaches in Prediction of Plant Virus Encoded miRNAs and Their Targeted Plant Genes.

This chapter presents a comprehensive approach to predict novel miRNAs encoded by plant viruses and identify their target plant genes, through integration of various ab initio computational approaches. The predictive process begins with the analysis of plant viral sequences using the VMir Analyzer software. VMir Viewer software is then used to extract primary hairpins from these sequences. To distinguish real miRNA precursors from pseudo miRNA precursors, MiPred web-based software is employed. Verified real pre-miRNA sequences with a minimum free energy of < -20 Kcal/mol, are further analyzed using the RNAshapes software. Validation of predictions involves comparing them with available Expressed Sequence Tags (ESTs) from the relevant plant using BlastN. Short sequences with lengths ranging from 19 to 25 nucleotides and exhibiting <5 mismatches are prioritized for miRNA prediction. The precise locations of these short sequences within pre-miRNA structures generated using RNAshapes are meticulously identified, with a focus on those situated on the 5' and 3' arms of the structures, indicating potential miRNAs. Sequences within the arms of pre-miRNA structures are used to predict target sites within the ESTs of the specific plant, facilitated by psRNA Target software, revealing genes with potential regulatory roles in the plant. To confirm the outcome of target prediction, results are individually submitted to the RNAhybrid web-based software. For practical demonstration, this approach is applied to analyze African cassava mosaic virus (ACMV) and East African cassava mosaic virus-Uganda (EACMV-UG) viruses, as well as the ESTs of Jatropha and cassava.

Gene Co-expression Network Analysis to Mine Genes Involved in Plant Secondary Metabolite Biosynthesis and Regulation.

Plants produce diverse specialized metabolites (SMs) that do not participate in plant growth and development but help them adapt to various environmental conditions. In addition to aiding in plant adaptation, different SMs serve as active ingredients for pharmaceutical and cosmetics products. However, despite their significant role in plant adaptation and industrial importance, the genes involved in the biosynthesis and regulation of many SMs remain largely unknown. This hinders deciphering the specific role of SMs in plant adaptation and limits their industrial utilization. Since many SMs pathway genes are expected to act in tight association with each other within a coexpression network, the network biology approach, such as weighted gene coexpression network analysis, could be used to identify the unknown genes. This chapter describes a workflow for constructing a gene coexpression network to identify genes that could be associated with the biosynthesis and regulation of SMs.

Characterization and identification of sources of rust resistance in Triticum militinae derivatives.

Triticum militinae (2n = 4X = 28, AtAtGG), belonging to the secondary gene pool of wheat, is known to carry resistance to many diseases. Though some disease resistance genes were reported from T. timopheevii, the closest wild relative of T. militinae, there are no reports from T. militinae. Twenty-one T. militinae Derivatives (TMD lines) developed at the Division of Genetics, IARI, New Delhi, were evaluated for leaf and stripe rusts at seedling and adult plant stages. Eight TMD lines (6-4, 6-5, 11-6, 12-4, 12-8, 12-12, 13-7 and 13-9) showed seedling resistance to both leaf and stripe rusts while six TMD lines (7-5, 7-6, 11-5, 13-1, 13-3 and 13-4) showed seedling resistance to leaf rust but adult plant resistance to stripe rust and three TMD lines (9-1, 9-2 and 15) showed seedling resistance to leaf rust but susceptibility to stripe rust. Three TMD lines (2-7, 2-8 and 6-1) with adult plant resistance to leaf and stripe rusts were found to carry the known gene Lr34/Yr18. Ten TMD lines (7-5, 7-6, 9-1, 9-2, 11-5, 11-6, 12-12, 12-4, 12-8, and 15) with seedling resistance to leaf rust, showing absence of known genes Lr18 and Lr50 with linked markers requires further confirmation by the test of allelism studies. As not a single stripe rust resistance gene has been reported from T. militinae or its close relative T. timpopheevii, all the 8 TMD lines (6-4, 6-5, 11-6,12-4, 12-8, 12-12, 13-7 and 13-9) identified of carrying seedling resistance to stripe rust and 3 TMD lines (13-1, 13-3 and 13-4) identified of carrying adult plant resistance to stripe rust are expected to carry unknown genes. Also, all the TMD lines were found to be cytologically stable and thus can be used in inheritance and mapping studies.

Systematic characterization of Gossypium GLN family genes reveals a potential function of GhGLN1.1a regulates nitrogen use efficiency in cotton.

The enzyme glutamine synthetase (GLN) is mainly responsible for the assimilation and reassimilation of nitrogen (N) in higher plants. Although the GLN gene has been identified in various plants, there is little information about the GLN family in cotton (Gossypium spp.). To elucidate the roles of GLN genes in cotton, we systematically investigated and characterized the GLN gene family across four cotton species (G. raimondii, G. arboreum, G. hirsutum, and G. barbadense). Our analysis encompassed analysis of members, gene structure, cis-element, intragenomic duplication, and exploration of collinear relationships. Gene duplication analysis indicated that segmental duplication was the primary driving force for the expansion of the GhGLN gene family. Transcriptomic and quantitative real-time reverse-transcription PCR (qRT-PCR) analyses indicated that the GhGLN1.1a gene is responsive to N induction treatment and several abiotic stresses. The results of virus-induced gene silencing revealed that the accumulation and N use efficiency (NUE) of cotton were affected by the inactivation of GhGLN1.1a. This study comprehensively analyzed the GhGLN genes in Gossypium spp., and provides a new perspective on the functional roles of GhGLN1.1a in regulating NUE in cotton.

Integrated genome-wide association and transcriptomic analysis to identify receptor kinase genes to stripe rust resistance in wheat germplasm from southwestern China.

Stripe rust of wheat, caused by Puccinia striiformis f. sp. tritici (Pst), is one of the most important diseases of wheat worldwide. Identification of new and elite Pst-resistance loci or genes has the potential to enhance overall resistance to this pathogen. Here, we conducted an integrated genome-wide association study (GWAS) and transcriptomic analysis to screen for loci associated with resistance to stripe rust in 335 accessions from Yunnan, including 311 landraces and 24 cultivars. Based on the environmental phenotype, we identified 113 protein kinases significantly associated with Pst resistance using mixed linear model (MLM) and generalized linear model (GLM) models. Transcriptomic analysis revealed that 52 of 113 protein kinases identified by GWAS were up and down regulated in response to Pst infection. Among these genes, a total of 15 receptor kinase genes were identified associated with Pst resistance. 11 candidate genes were newly discovered in Yunnan wheat germplasm. Our results revealed that resistance alleles to stripe rust were accumulated in Yunnan wheat germplasm, implying direct or indirect selection for improving stripe rust resistance in elite wheat breeding programs.

Functional analysis of a wheat class III peroxidase gene, TaPer12-3A, in seed dormancy and germination.

BACKGROUND: Class III peroxidases (PODs) perform crucial functions in various developmental processes and responses to biotic and abiotic stresses. However, their roles in wheat seed dormancy (SD) and germination remain elusive. RESULTS: Here, we identified a wheat class III POD gene, named TaPer12-3A, based on transcriptome data and expression analysis. TaPer12-3A showed decreasing and increasing expression trends with SD acquisition and release, respectively. It was highly expressed in wheat seeds and localized in the endoplasmic reticulum and cytoplasm. Germination tests were performed using the transgenic Arabidopsis and rice lines as well as wheat mutant mutagenized with ethyl methane sulfonate (EMS) in Jing 411 (J411) background. These results indicated that TaPer12-3A negatively regulated SD and positively mediated germination. Further studies showed that TaPer12-3A maintained H2O2 homeostasis by scavenging excess H2O2 and participated in the biosynthesis and catabolism pathways of gibberellic acid and abscisic acid to regulate SD and germination. CONCLUSION: These findings not only provide new insights for future functional analysis of TaPer12-3A in regulating wheat SD and germination but also provide a target gene for breeding wheat varieties with high pre-harvest sprouting resistance by gene editing technology.

Identification of QTNs, QTN-by-environment interactions, and their candidate genes for salt tolerance related traits in soybean.

BACKGROUND: Salt stress significantly reduces soybean yield. To improve salt tolerance in soybean, it is important to mine the genes associated with salt tolerance traits. RESULTS: Salt tolerance traits of 286 soybean accessions were measured four times between 2009 and 2015. The results were associated with 740,754 single nucleotide polymorphisms (SNPs) to identify quantitative trait nucleotides (QTNs) and QTN-by-environment interactions (QEIs) using three-variance-component multi-locus random-SNP-effect mixed linear model (3VmrMLM). As a result, eight salt tolerance genes (GmCHX1, GsPRX9, Gm5PTase8, GmWRKY, GmCHX20a, GmNHX1, GmSK1, and GmLEA2-1) near 179 significant and 79 suggested QTNs and two salt tolerance genes (GmWRKY49 and GmSK1) near 45 significant and 14 suggested QEIs were associated with salt tolerance index traits in previous studies. Six candidate genes and three gene-by-environment interactions (GEIs) were predicted to be associated with these index traits. Analysis of four salt tolerance related traits under control and salt treatments revealed six genes associated with salt tolerance (GmHDA13, GmPHO1, GmERF5, GmNAC06, GmbZIP132, and GmHsp90s) around 166 QEIs were verified in previous studies. Five candidate GEIs were confirmed to be associated with salt stress by at least one haplotype analysis. The elite molecular modules of seven candidate genes with selection signs were extracted from wild soybean, and these genes could be applied to soybean molecular breeding. Two of these genes, Glyma06g04840 and Glyma07g18150, were confirmed by qRT-PCR and are expected to be key players in responding to salt stress. CONCLUSIONS: Around the QTNs and QEIs identified in this study, 16 known genes, 6 candidate genes, and 8 candidate GEIs were found to be associated with soybean salt tolerance, of which Glyma07g18150 was further confirmed by qRT-PCR.

A combination of QTL mapping and genome-wide association study revealed the key gene for husk number in maize.

KEY MESSAGE: Two key genes Zm00001d021232 and Zm00001d048138 were identified by QTL mapping and GWAS. Additionally, they were verified to be significantly associated with maize husk number (HN) using gene-based association study. As a by-product of maize production, maize husk is an important industrial raw material. Husk layer number (HN) is an important trait that affects the yield of maize husk. However, the genetic mechanism underlying HN remains unclear. Herein, a total of 13 quantitative trait loci (QTL) controlling HN were identified in an IBM Syn 10 DH population across different locations. Among these, three QTL were individually repeatedly detected in at least two environments. Meanwhile, 26 unique single nucleotide polymorphisms (SNPs) were detected to be significantly (p < 2.15 × 10-6) associated with HN in an association pool. Of these SNPs, three were simultaneously detected across multiple environments or environments and best linear unbiased prediction (BLUP). We focused on these environment-stable and population-common genetic loci for excavating the candidate genes responsible for maize HN. Finally, 173 initial candidate genes were identified, of which 22 were involved in both multicellular organism development and single-multicellular organism process and thus confirmed as the candidate genes for HN. Gene-based association analyses revealed that the variants in four genes were significantly (p < 0.01/N) correlated with HN, of which Zm00001d021232 and Zm00001d048138 were highly expressed in husks and early developing ears among different maize tissues. Our study contributes to the understanding of genetic and molecular mechanisms of maize husk yield and industrial development in the future.

A wild rice CSSL population facilitated identification of salt tolerance genes and rice germplasm innovation.

Salt stress is one of the major factors that limits rice production. Therefore, identification of salt-tolerant alleles from wild rice is important for rice breeding. In this study, we constructed a set of chromosome segment substitution lines (CSSLs) using wild rice as the donor parent and cultivated rice Nipponbare (Nip) as the recurrent parent. Salt tolerance germinability (STG) was evaluated, and its association with genotypes was determined using this CSSL population. We identified 17 QTLs related to STG. By integrating the transcriptome and genome data, four candidate genes were identified, including the previously reported AGO2 and WRKY53. Compared with Nip, wild rice AGO2 has a structure variation in its promoter region and the expression levels were upregulated under salt treatments; wild rice WRKY53 also has natural variation in its promoter region, and the expression levels were downregulated under salt treatments. Wild rice AGO2 and WRKY53 alleles have combined effects for improving salt tolerance at the germination stage. One CSSL line, CSSL118 that harbors these two alleles was selected. Compared with the background parent Nip, CSSL118 showed comprehensive salt tolerance and higher yield, with improved transcript levels of reactive oxygen species scavenging genes. Our results provided promising genes and germplasm resources for future rice salt tolerance breeding.

Genome-wide characterization of DNA methyltransferase family genes implies GhDMT6 improving tolerance of salt and drought on cotton.

BACKGROUND: DNA methylation is an important epigenetic mode of genomic DNA modification and plays a vital role in maintaining epigenetic content and regulating gene expression. Cytosine-5 DNA methyltransferase (C5-MTase) are the key enzymes in the process of DNA methylation. However, there is no systematic analysis of the C5-MTase in cotton so far, and the function of DNMT2 genes has not been studied. METHODS: In this study, the whole genome of cotton C5-MTase coding genes was identified and analyzed using a bioinformatics method based on information from the cotton genome, and the function of GhDMT6 was further validated by VIGS experiments and subcellular localization analysis. RESULTS: 33 C5-MTases were identified from three cotton genomes, and were divided into four subfamilies by systematic evolutionary analysis. After the protein domain alignment of C5-MTases in cotton, 6 highly conserved motifs were found in the C-terminus of 33 proteins involved in methylation modification, which indicated that C5-MTases had a basic catalytic methylation function. These proteins were divided into four classes based on the N-terminal difference, of which DNMT2 lacks the N-terminal regulatory domain. The expression of C5-MTases in different parts of cotton was different under different stress treatments, which indicated the functional diversity of cotton C5-MTase gene family. Among the C5-MTases, the GhDMT6 had a obvious up-regulated expression. After silencing GhDMT6 with VIGS, the phenotype of cotton seedlings under different stress treatments showed a significant difference. Compared with cotton seedlings that did not silence GhDMT6, cotton seedlings silencing GhDMT6 showed significant stress resistance. CONCLUSION: The results show that C5-MTases plays an important role in cotton stress response, which is beneficial to further explore the function of DNMT2 subfamily genes.

Information-incorporated gene network construction with FDR control.

MOTIVATION: Large-scale gene expression studies allow gene network construction to uncover associations among genes. To study direct associations among genes, partial correlation-based networks are preferred over marginal correlations. However, FDR control for partial correlation-based network construction is not well-studied. In addition, currently available partial correlation-based methods cannot take existing biological knowledge to help network construction while controlling FDR. RESULTS: In this paper, we propose a method called Partial Correlation Graph with Information Incorporation (PCGII). PCGII estimates partial correlations between each pair of genes by regularized node-wise regression that can incorporate prior knowledge while controlling the effects of all other genes. It handles high-dimensional data where the number of genes can be much larger than the sample size and controls FDR at the same time. We compare PCGII with several existing approaches through extensive simulation studies and demonstrate that PCGII has better FDR control and higher power. We apply PCGII to a plant gene expression dataset where it recovers confirmed regulatory relationships and a hub node, as well as several direct associations that shed light on potential functional relationships in the system. We also introduce a method to supplement observed data with a pseudogene to apply PCGII when no prior information is available, which also allows checking FDR control and power for real data analysis. AVAILABILITY AND IMPLEMENTATION: R package is freely available for download at https://cran.r-project.org/package=PCGII.